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Abstract We estimated the potential habitat for an introduced population of Hestina assimilis assimilis (Linnaeus) 
in Japan, based on climates and the potential habitat of three Celtis plants that are hosts or potential hosts. We showed 
that many areas in Japan had a suitable climate for H. a. assimilis. This indicates that 万 a. assimilis could live in 
various areas of Japan where Celtis plants are found. Celtis sinensis Pers. and Celtis boninensis Koidz. preferred 
an artificial environment, but C. sinensis had a wide range of suitable habitat. Our results indicate that the potential 
habitat of H. a. assimilis in Japan is much wider than the present distribution area of the introduced population. It is 
necessary to perform further impact assessment on the native ecosystem and to establish appropriate management of H. 
a. assimilis in Japan. 
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は じ め に 


生物 多様 性 を 劣化 させ る 要因 と し て , 外来 生物 の 分 布 拡大 
は 世界 規模 で の 問題 に な っ て いる . チョ ウ 類 に お いて も 本 
来生 息 し て いな か っ た 地域 へ 侵入 , 定着 する 機会 が 増え て 
いる . この よう な 人 外来 生物 の 侵入 経路 と し て は , 寄主 植物 
で ある 農作物 や 観賞 用 植物 , 街路 樹 な どの 栽培 や 植 栽 を 通 
じ て 広 が る よう な 非 意図 的 な 経路 や , 愛好 家 に よる 放 蝶 と 
いっ た 意図 的 な 経路 な どき さま ざま で ある . いずれ に お いて 


因 か ら 生 物 の 生息 適地 を 予測 する こと が 可能 で ある . し か 
L, 分 布 拡大 し て いる 種 を 解析 し よう と する 場合 , 分 布 情 
報 が 十分 に 蓄積 さん て いな いこ と も 多く , 適切 な 解析 が 困 
難 な 場 合 も ある . 侵入 初期 の 段階 で 対策 を 考え て いく た め 
に は , 情報 の 蓄積 を 待た ず に 解析 が 実施 で きる こと が 望ま 
NS. この よう な 状況 の 場合 , 侵入 ・ 分 布 拡大 後に 鞭 積 さ 
れ た 情報 を 用 いる の で は な く , 原産 地 の 分 布 情報 を 気候 な 
どの 環境 要因 か ら 予 測 す る モデ ル を 構築 し その モデ ル に 
侵入 先 な どの 当該 地域 の 環境 要因 を 与え る こと で , その 地 


も , MUD EAL, 在 来 種 と の 競合 や 交雑 な どの 危険 性 
を は らん で お り , 侵入 地域 の 生態 系 に 大 き な 影 響 を 与え る 
可能 性 が ある RA + AEF, 2003 : 矢 後 , 2012). 既存 の 生 
能 系 へ の 影響 を 最小 限 に する に は , DHEA ZIRHO 
適切 な 管理 を 考え て いく 必要 が あり , その た め に は 将来 の 
分 布 を 見 積もる こと が 重要 で ある . 

生物 の 将来 の 分 布 を 知る た め に は , 種 の 分 布 モ デル を 
用 いて 潜在 的 な 生息 適地 を 解析 する こと が 有効 で ある 
(Franklin, 2009), 種 の 分 布 モ デル は 生物 の 分 布 と 気候 や 土 
地 利 用 な どの 要因 と の 関係 を 説明 する も の で , それ ら の 要 
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域 の 潜在 的 な 生息 適地 を 抽出 する 手法 を 用 いる こと が 可能 
で ある (例え ば , オオ クチ バス ・ コ クチ バス : Iguchi etal., 
2004, ヒト スジ シマ カ :Medley, 2010). この 手法 に よっ て . 
侵入 前 や 侵入 初期 の 段階 に お ける 外来 種 の リス ク 管 理 に 対 
し て 重要 な 知見 を 提供 する こと が で きる . た だ し , 侵入 地 
域 の 中 に 原産 地 に 存在 し な い 環 境 (つま り , その 種 に と っ 
て 未経験 の 環境 ) が 含ま れる 場合 、 この よう な アプ ロー チ 
で は 推定 結果 が 信頼 で き な い こと が 指摘 され て いる の で 注 
意 が 必要 で ある (Too etal., 2007). 
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日 本 に お いて は , 近年 , アカ ボ シ ゴ マダ ラ の 分 布 拡大 が 
これ まで 分 布 し て いな か っ た 関東 地方 を 中 心 に 報告 され 
て いる (岩野 , 2005 : 岩野 ・ 菅 井 , 2007 : 月 刊 む し 編集 部 , 
2010 : 柄 澤 ほ か , 2012). 日 本 に お いて は , 在 来 の アカ ボ シ 
ゴマ ダラ と し て 奄美 諸島 に 生息 する 日 本 固有 亜 種 Hesrina 
assimilis shirakii Shir6zu が 知ら れ て いる が , 関東 地方 で 見 
られ る アカ ボ シ ゴ マダ ラ は その 形態 的 な 特徴 か ら 大 陸 産 
の 名 義 タ イプ 亜 種 万 es assimilis assimilis (Linnaeus) で あ 
る と され る (浅野 , 2002a, 2002b : 小野 寺 , 2002 : 中 村 ほ か , 
2003) (AA CALE, H. a shirakii を 在 来 アカ ボ シ ゴ マ 
ダラ , H. a. assimilis を 外来 アカ ボ シ ゴ マダ ラ と 表現 する ). 
この 外来 アカ ボ シ ゴ マダ ラ の 侵入 し た 経緯 に つい て は , H 
蝶 に よる も の で ある 可能 性 が 高い と 考え られ て いる (HH: 
菅井 , 2005 ; 松井 , 2010 ; 高桑 ,2012). 外来 アカ ボ シ ゴ マダ 
ラ は , 効 虫 期 に エノキ 属 植物 を 餌 と する た め , 寄主 植物 が 
同じ で ある オオ ムラ サキ や ゴマ ダラ チョ ウ と の 資源 競合 や 
ニッ チ ナチ 競 合 が 危 選 きれ て いる (松井 , 2010; 矢 後 , 2012). 
また , 奄美 諸島 に 外来 アカ ボ シ ゴ マダ ラ が 侵入 し た 場合 , 
在 来 ア カボ シゴ マダ ラ と の 交配 に よる 遺伝 的 汚染 の 危険 性 
が 考え られ る . さら に , 在 来 の 近 緑 種 ゴ マダ ラチ ョ ウ と 
の 野外 環境 で の 交雑 (石垣 . 2000) や オオ ムラ サキ と の 交 
尾 (大 石 , 2011) も 確認 され て お り , 他 種 ま で 含め た 遺伝 
的 失 乱 や 繁殖 阻害 ・ 和 干渉 も 懸念 され て いる (高桑 , 2012) . 
この よう な 懸念 が ある た め , 外来 アカ ボ シ ゴ マダ ラ は 人 外来 
生物 法 に よっ て 要 注 意外 来生 物 に 指定 され て いる (環境 省 , 
2014b). 


外来 アカ ボ シ ゴ マダ ラ は 分 布 拡大 途中 に あり , 今後 も 分 布 
を 広げ て いく も の と 考え られ る . また , 本 種 は 何ら か の 
事情 で 神奈 川 県 か ら 侵入 ・ 分 布 拡大 し た と 推測 され て いる 
が , 国内 に お ける 他 の 地方 に も 侵入 ・ 定 着 す る 可能 性 も あ 
る . 実際 , 関東 地方 か ら 離 れ た 京都 や 名 古屋 で も 外来 アカ 
ボ シ ゴ マダ ラ が 確認 され て いる (間野 ・ 岩 本 , 2011 中 橋 , 
2011 ; 渡邊 . 2013a, 2013b). 一 般 的 に は , 生物 の 分 布 は 気 
候 に よる 制約 を 受け る た め , 日 本 に お ける 外来 アカ ボ シ ゴ 
マダ ラ に お いて も 気候 に よる 制約 を 考慮 し た 上 で その 洪 在 
的 な 生息 適地 を 把握 する 必要 が ある . 日 本 に お ける 外来 ア 
カボ シゴ マダ ラム の 分 布 情報 の み で は その 収集 地域 が 限ら れ 
て いる た め , 解析 に 用 いる 気候 の 幅 が 狭く , 広域 へ の 外 挿 
に 耐え られ る よう な 信頼 に 足る 解析 結果 を 得 ら れ な い 可 能 
性 が ある (Loo et al., 2007). し か し , 本 来 の 分 布地 域 で あ 
る 大 陸 に お ける 標本 情報 を 活用 すれ ば , 採集 地域 が 広く 十 
分 な 気候 の 幅 を 踏ま えて 解析 で きる の で , 気候 の 視点 か ら 
見 た 外来 アカ ボ シ ゴ マダ ラ の 生息 適地 を 解析 する こと が 可 
能 で ある . 一 方 , チョ ウ 類 の 分 布 に は 寄主 植物 の 存在 が 必 
要 不 可 欠 で ある . 外来 アカ ボ シ ゴ マダ ラ が 生息 可能 な 場所 
は , 気候 が 適し て いる と いう こと 以外 に , 寄主 植物 で ある 
エノキ 属 植物 が 分 布 し て いる こと が 必要 で ある . そこ で , 
エノキ 属 植物 の 分 布 情報 を 用 いて , エノキ 属 植物 の 生息 適 
地 解析 も お こ な う . これ に より , 寄主 植物 の 観点 か ら 見 た 
外来 アカ ボ シ ゴ マダ ラ の 潜在 的 な 生息 適地 を 明らか に する 
こと が で きる . その 上 で 気候 に よる 外来 アカ ボン ゴマ ダラ 
の 生息 適地 と 寄主 植物 の 生息 適地 を 合わ せ て 考察 する こと 


で , 本 種 の 潜在 的 な 生息 適地 の 把握 が 可能 と 考え られ る . 


本 研究 で は , チョ ウ の 標本 情報 と 寄主 植物 の 分 布 情報 を 用 
いて , 日 本 に お ける 外来 アカ ボ シ ゴ マダ ラ の 潜在 的 な 生息 
適地 を 推定 する こと を 目的 と する . 本 研究 に お ける 成果 は 
さま ざま な 外来 種 に お ける 管理 指針 の 策定 へ の 寄与 が 期待 
され る は ほか , 他 の 外来 器 に お ける 解析 へ の 応用 も 期待 され 
る . 


材料 お よび 方 法 

1. 気候 要因 に 基づく 外来 アカ ボ シ ゴ マダ ラ の 浴 在 的 生息 適 
地 解析 

まず , 標本 の 位置 情報 が 掲載 され て いる 文献 (Lang, 
2012 : 増井 ・ 玉 井 , 2013) を 用 いて , 中 国 お よび 韓国 に お 
ける 外来 アカ ボ シ ゴ マダ ラ の 分 布地 点 の 緯度 経度 を 入手 し 
Zz, これ ら の 地域 で は 万 a. assimilis は 普通 種 で あり , 通年 
生息 し て いる と 考え られ る (Park and Kim, 1997: 増 井 ・ 玉 井 , 
2013). Lang (2012) に お いて は , 1930 年 代 か ら 2000 年 代 
に 得 ら れ た a. gssz が s と され て いる 標本 情報 か ら 34 地 点 
の 位置 情報 を 利用 し , 他 の 亜 種 に され て いる 標本 情報 は 利 
用 し な か っ た . 増井 ・ 玉 井 (2013) に つい て は , Hestina 
assimilis inexpecta Masui & Tamai を 除い て , 記さ れ て いた 15 
地点 の 位置 情報 を 使用 し た . これ ら を 合わ せ て 49 地 点 の 位 
置 情報 を 入手 し た (Mla). 


外来 アカ ボ シ ゴ マダ ラ の 分 布 を 説明 する た め の 気 候 要 因 と 
LT, 中 国 ・ 韓 国 お よび 日 本 に お ける 最 寒 月 の 最低 気温 ( 以 
T, 最低 気温 ) と 冬季 (12 月 ~ 3 月 ) の 最低 月 平均 相対 湿 
E (以下 , 冬季 最低 湿度 ) を , 位置 情報 を 持つ デー タ を 統 
合 ・ 管 理 す る シス テム で ある GTIS (Geographic Information 
System : 地理 情報 シス テム ) を 用 いて 用 意 し た . 冬季 の 
最低 気温 や 湿度 (すなわち 乾燥 の 程度 ) は , 外来 アカ ボ 
シゴ マダ ラ の 死亡 率 に 寄与 する 重要 な 要因 で ある と 考え 
た . 最低 気温 は , Bioclim (Hijimans et al., 2014) に 用 意 き 
れ て いる 紳 度 経度 解像度 10 分 の BIO6 を 用 いた ( 図 tb). Z 
れ は jjimans et al. (2005) に 基づい て 計算 され た 1950 年 か 
ら 2000 年 に お ける 気候 ダー タ か ら 作 成 さ れ た デー タ セ ッ ト 
で ある . 使用 し た 標本 情報 の 年 代 と 合わ せる た め に , Ait 
究 で は この デー タ を 使用 し た . 冬季 最低 湿度 に は , 全球 
規模 の 年 平均 相対 湿度 (解像度 5 分 ) の デー タ (New etal, 
1999) を 使用 し た . KEL, この まま で は 冬季 の 湿度 が 分 
か ら な いこ と か ら , NASA (2007) に よっ て 作成 され た 解 
像 度 の 粗い (1 度 ) 月 ご と の 平均 相対 湿度 分 布 デ ー タ を 用 
いて 年 平均 湿度 に 対す る 冬季 最低 月 平均 気温 の 比 を 算出 
L, それ を 逆 距離 加重 法 (Inverse Distance Weighting) に よっ 
て 内 挿 補間 し , New etal. (1999) の 年 平均 相対 湿度 に 乗 
じ る こと で 冬季 最低 湿度 を 算出 し た (Mic). な お , 最低 
気温 と 冬季 最低 湿度 の 解像度 は 10 分 に 統一 し た . 


桂 本 情報 に 基づく 分 布地 点 情 報 は , 種 の 分 布 解析 に よく 用 
いら れる が , ある 場所 に 「 い な い 」 COIR Ted 
ちと 不在 地点 の 情報 が ほとん どの 場合 存在 し な いと いう 欠点 
が ある . その た め , 分 布 解析 に は 最大 エン トロ ビー 法 に 基 
づく 統計 モデ リン グ 手 法 で あり , 在 地点 の みか ら で も 比較 
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a) 標本 情報 に 基づく アカ ボ シ ゴ マダ ラ の 分 布 











解析 対象 範囲 











図 1. 文献 に お ける 標本 情報 か ら 抽 出し た アカ ボ シ ゴ マダ ラ の 位 
置 情 報 . (。) お よび 最 寒 月 の 最低 気温 (b), 4% (12 月 一 3 月 ) 
の 最低 月 平均 相対 湿度 (c) の 空間 分 布 . 


的 頑健 な 解析 が 可能 な Maxent (Phillips et al., 2006 : Elith et 
al., 2011) を 採用 し た . 解析 ソフ トウ ェ ア に は , 同 手法 に 
特 化し た Maxent バ ー ジ ョ ン 3.3.3k (Dudík et al., 2011) を 用 
いた . モデ ル 構 築 の 対象 範囲 は , 中 国 お よび 韓国 と し た 
(Mla). モデ ル の オー バー フィ ッ テ ィ ング を 避け る た め に , 
モデ ル の 形 は 1 次 項 お よび 2 次 項 を 使用 し た . モデ ル の 当 
て は まり は , モデ ル 全 体 の 当て は まり を 示す こと が で きる 
ROC (Receiver Operating Characteristics) 曲線 の AUC (Area 
Under Curve) の 値 で 評価 し た (Fielding and Bell, 1997). 
中 国 お よび 韓国 を 対象 と し て Maxent に よっ て 構築 し た モデ 
ル を 日 本 に 外 挿す る こと で , 日 本 に お ける 気候 要因 に 基 づ 
く 外 来 ア カボ シゴ マダ ラ の 潜在 的 生息 適地 を 算出 し た . 一 
連 の 解析 手順 を ., ブー トス トラ ッ プ 法 で リサ ンプ リン グ を 
200 回 繰り 返す こと に よっ て , モデ ル の 信頼 性 を 評価 し た . 


2. 寄主 植物 の 潜在 的 生息 適地 解析 


外来 アカ ボ シ ゴ マダ ラ の 寄主 植物 は エノキ 属 の 樹木 で ある 
(日 本 チョ ウ 類 保全 協会 . 2012 : 国立 環境 研究 所 , 2013). 





その た め , エノキ 属 植物 の 分 布 情報 を 得る た め に , 環境 省 
に よる 第 6 回 ・7 回 自然 環境 保全 基礎 調査 に お いて 広域 的 に 
実施 され た 植生 調査 結果 を 用 いた (環境 省 , 2014a). この 
調査 で は 日 本 を 8 つの 区 域 に 分 割 し . それ ぞ れ の 区 域 で 調 
査 が お こ な わ れ て お り , 本 研究 で は 2000 年 か ら 2010 年 まで 
に 得 ら れ た 14.688 地 点 の 調査 結果 を 利用 し た . た だ し , 後 
述 す る 環境 要因 が 揃わ な か っ た 地点 は 解析 か ら 除 いた た 
め , 実際 に 使用 し た 地点 は 14.428 地 点 で ある . この デー 
タ か ら , 出現 地点 数 が 多かっ た 3 種 (エノキ Celtis sinensis 
Pers, エゾ エノキ Celtis jessoensis Koidz, クリ ノ ハ エノキ 
Celtis boninensis Koidz.) を 解析 に 用 いた . それ ぞ れ の 出現 
地点 数 は , 863, 113, 31 地 点 で あっ た . 


それ ぞ れ の 種 に つい て , 植生 調査 地点 ご と に 出現 の 有無 
(10) デー タ に 変換 し . その 出現 の 有無 を 環境 要因 で 説明 
する 分 市 も デル を 構 条 し た 、 ただし, エノキ と クワ シメ ハエ 
ノ キ は 北海 道 に は 自生 し て いな いこ と か ら , これ ら 2 種 に 
つい て は 北海 道 に お ける 植生 調査 地点 を 解析 対象 か ら 除 
外し た . 環境 要因 に は , 3 次 メッ シュ で 用 意 し た 気候 と 地 
É, 土地 利用 を 用 いた . 気候 要因 に は , 温 量化 指数 (Kira, 
1991), 最 寒 月 最低 気温 (Ohsawa, 1993), 夏季 ($ 月 一 9 月 ) 
隆 水 量 お よび 冬季 (12 月 一 3 月 ) 降水 量 を 用 いた . 上 記 の 
気候 要因 の 算出 に は , 3 次 メッ シュ の 値 が 用 意 さ れ て いる 
メッ シュ 気候 値 2000 を 用 いた . 地形 要因 に は , 地形 的 湿潤 
数 お よび 地形 区 分 を 用 いた . 地形 的 湿潤 指数 は 斜面 位 軒 
に よっ て 生じ る 土壌 水分 傾 度 を 数 量 的 に あら わす 指数 で あ 
る (Moore etal., 1993 ; Sorensen et al., 2006). 地形 区 分 は . 
若松 ほか (2005) に お ける 3 次 メッ シュ ご と の 地形 区 分 を 
低地 . 丘陵 地 お よび 山地 の 3 つ に 再 分 類 し た も の を 用 いて ., 
丘陵 地 1) と それ 以外 (0), 山地 (1) と それ 以外 (0) 
と いう 2 つの ダミ ー 変 数 を 作成 し . それ ぞ れ 地形 の 変数 と 
し た . 土地 利用 要因 に は 人 間 活 動 と 関連 の 深い , 農地 率 お 
よび 都市 率 を 用 い , さら に 里山 地域 な ど に お ける 土地 利用 
の 多様 性 も 重要 で ある 可能 性 を 考慮 し . 景観 の 不 均一 性 も 
算出 し た . これ ら 土 地 利 用 と 景観 不 均一 性 の 値 は , 周辺 2 
メッ シュ まで を 含め て 値 を 算出 し た . ベー ス と な る 土地 利 
用 デー タ は , 環境 省 の 第 2-$ 回 自然 環境 保全 基礎 調査 に よ 
る 植生 図 を 再 分 類 し て 作成 し た (小川 ほか . 2013 : Akasaka 
et al., in press). 景観 不 均一 性 は , Kadoya and Washitani(2011) 
に よる Satoyama Index を 補正 し て 算出 し た . Kadoya and 
Washitani (2011) で は 都市 的 土地 利用 が 優 占 する 場所 を 解 
析 対 象 か ら 除 外し て いた が . 本 人 研究 で は それ も 含め て 値 を 
算出 し た . 


分 布 モ デル に は , ロジ ステ ィ ッ ク 回 帰 分 析 を 用 いた . これ 
に より , 分 布 確率 を 予測 する こと が で きる . この 解析 に は 
R バ ー ジ ョ ン 3.0.2 (R Core Team, 2013) に お ける stats パ ッ 
ケー ジ の glm 関 数 を 用 いた . 種 ご と に , すべ て の 変数 の 組 
み 合 わせ で モデ ル の 当て は まり の 指標 で ある AIC (Akaike 
Information Criterion) を 算出 し , 最も AIC が 小さ か っ た 組 
み 合 わせ を 最適 モデ ル と し た . モデ ル の 精度 の 評価 に は 
AUC を 用 いた . また ., 分 布 の 適地 と 不適 地 を 分 ける た め に , 
モデ ル の 感度 と 特異 度 の 合計 値 が 最大 と な る 予測 値 を 隊 値 
と し て 算出 し た . 得 ら れ た 最適 モデ ル を 日 本 全体 に 当て は 
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め る こと で , 面 的 な 分 布 確率 に 基づく 人 息 適地 推定 を お こ 
Rok. な お , エノキ と クワ リノ ハ エノキ は モデ ル 作 成 時 に 
は 北海 道 を 対象 外 と し た が , 外 挿 の 際 に は 北海 道 も 含め た . 


結 R 


1. 気候 要因 に 基づく 外来 アカ ボ シ ゴ マダ ラ の 凌 在 的 生息 敵 
地 

Maxent に よっ て 得 ら れ た モデ ル の AUC は 0.84 よ 0.02 (平均 
値 + 標 準 誤差 ) で あっ た こと か ら , 良い 精度 を 持っ た モデ 
ル を 構築 する こと が で きた . モデ ル に よる 気候 要因 の 反応 
曲線 に よれ ば , 中 国 お よび 韓国 に お ける 外来 アカ ボ シ ゴ マ 
ダラ の 生息 適地 は , 最低 気温 お よび 冬季 最低 湿度 が 高い 地 
域 で あっ た ( 図 2). 地図 上 で みる と , 大 陸 の 南東 部 で 生息 
適合 度 が 高い こと を 確認 で きる ( 図 3). 得 ら れ た モデ ル を 
日 本 に 外 挿し た 結果 , 日 本 の ほとん どの 地域 で 生息 適合 度 
が 高い 傾向 が 見 られ た (3). 


2. 寄主 植物 の 潜在 的 生息 適地 


エノキ 属 3 種 に つい て ロジ ステ ィ ッ ク 回 帰 分 析 を お こない . 
それ ぞ れ の 最適 モデ ル の AUC を 調べ た 結果 . エノキ で は 
0.72, エゾ エノキ で は 0.67,. クワ ノ ハ エノキ で は 0.98 で あっ 
た . クワ ノ ハ エノキ は 精度 の 良い モデ ル が , エノキ で は ま 
ず ま ず の 精度 を 持っ た モデ ル が 得 ら れ た 一 方 で , エゾ エ ノ 
キ に つい て は あま り 精 度 の よい モデ ル は 得 ら れ な か っ た . 
得 ら れ た モデ ル を 見 て みる と , すべ て の 種 で 最低 気温 が 選 
択 さ れ , いずれ の 種 に お いて も 最低 気温 の 高い 地域 を 好 ん 
で いた . 特に , クワ ノ ハ エノキ は 最低 気温 が 最も 重要 な 変 
数 で あっ た . また , エノキ お よび クワ リノ ハエ ノ キ は 都市 率 
や 農地 率 と いっ た 人 為 的 な 環境 を 示す 要因 が 重要 な 変数 と 
し て 選択 され , それ ら の 環境 を 好 ん で いる と いう 結果 が 得 
られ た ( 表 1). 一 方 で エゾ エノキ に つい て は , ABN 
影響 に 関す る 変数 は 選択 され な か っ た . また , エゾ エノキ 
の 分 布 確率 は 他 の エノキ 類 の 分 布 確率 より も 著しく 低かっ 
た ( 図 4). 


考 g 


Maxent に よる 気候 要因 の 解析 か ら , 日 本 の 多く の 地域 に お 
いて , 気候 的 に は 外来 アカ ボ シ ゴ マダ ラ の 生息 が 可能 で あ 
る こと が 示唆 され た ( 図 3). 原産 地 の 気 候 を 用 いた 解析 を 
他 の 地域 に その まま 外 挿 する こと は , 未経験 の 気候 へ の 適 
応 の 評価 の 観点 に お いて 問題 生じ る こと も ある が (Loo 
et al., 2007), 本 研究 の 場合 , 日 本 の 気候 の 範囲 は 中 国 お よ 
び 韓 国 の 気候 の 範囲 に 含ま れる た め 大 き な 問 題 は 生じ な い 
と 考え られ る . 外来 アカ ボ シ ゴ マダ ラ は 最低 気温 の 高い 地 
域 に 生息 する 傾向 に あっ た が ( 図 2), 気温 の 低い 中 国 北部 
な ど で 採 集 さ れ た 標本 が 存在 し て お り (Lang, 2012), 気温 
か ら の 視点 と し て 日 本 の 広域 に 生息 可能 と いう 結果 は 妥当 
で ある と 考え られ る . また , オオ ムラ サキ で は 越冬 成功 の 
た め の 重 要 な 要因 の ひと つと し て 相葉 に 対す る 水分 供給 量 
の 多 さ が 指摘 きれ て お り (Kobayashi and Inaizumi, 2003), 
冬季 湿度 の 高い 地域 に お いて 外来 アカ ボ シ ゴ マダ ラ の 生 
息 適合 度 が 高い 理由 も これ と 同様 に 説明 で きる と 考え られ 
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図 2. Maxent に よっ て 構築 し た アカ ボ シ ゴ マダ ラ の 分 布 モ デル に 
お ける , 最 案 月 の 最低 気温 (a) お よび 冬季 (12 月 3 月 ) の 
最低 月 平均 相対 湿度 (b) の 反応 曲線 . 実線 は ブー トス トラ ッ 
TE (200 回 リサ ンプ リン グ ) に よる 中 央 値 , 点線 は 95% 信 頼 
区 間 を 示す . 


る . 気温 と 同様 に 原産 地 北 限 周辺 の 冬季 湿度 の 低い 地域 へ 
の 分 布 が 見 られ る こと か ら , 日 本 の 冬季 湿度 は 越冬 の 妨げ 
に は な ら な いと 考え られ る . 本 研究 の 結果 は 気候 要因 に よ 
る 外来 アカ ボ シ ゴ マダ ラ の 溶 在 分 布 適地 を お お ま か に 示せ 
た と 考え られ る . た だ し , 先行 研究 で は 休眠 誘起 や 発育 と 
いっ た 生活 史 の 観点 か ら , 外来 アカ ボン シ ゴマ ダラ は 日 本 の 
温帯 域 に 適応 し きれ て いな いと いう 指摘 が ある (加藤 ・ 宮 
内 , 2008). 現状 で は 関東 地方 に お ける 外来 アカ ボ シ ゴ マダ 
ラ は 分 布 拡大 を 続け て いる こと か ら 日 本 の 温帯 域 の 気候 で 
も ある 程度 集団 を 維持 で きる と 考え られ る が , 野外 に お け 
る 生存 や 成長 と 気候 と の 関係 や それ ら を 考慮 し た 浴 在 生息 
適地 の 把握 に つい て は 今後 の 課題 で あろ う . 


寄主 植物 で ある エノキ の 生息 適地 は 低地 や 平野 部 な どの 人 
為 的 な 影響 の 大 きい 場所 を 中 心 に 日 本 に 広く 分 布 し て お 
り , まだ 外来 アカ ボ シ ゴ マダ ラ の 定着 が 確認 され て いな い 
地域 に も 生息 適地 が 存在 し て いる (R3). エノキ は 陽 樹 的 
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図 3. Maxent に よっ て 構築 し た アカ ボ シ ゴ マダ ラ の 分 布 モ デル に 基づく アカ ボ シ ゴ マダ ラ の 生息 適合 度 . 
示さ れ て いる 値 は ブー トス トラ ッ プ 法 (200 回 リサ ンプ リン グ ) に よる 中 央 値 で ある . 


な 特徴 を 持つ 先駆 的 な 樹種 で あり , 措 乱 地 や 林 縁 部 , 都 
市 景観 な ど びに よく 出現 する こと が 知ら れ て いる ( 達 ・ 大 
IR, 1992 : 小林 ほか . 2004). 本 研究 の 結果 も その よう な 傾 
向 を 反映 する も の で あっ た . また , 北海 道 に は エノキ は 自 
生 し て いな い が , 植 栽 し た 場合 に 定着 で きる と 考え られ る 
地域 が 低地 の 都市 部 を 中 心 に 存在 し て いた . 暖か い 地 域 を 
好む な クワ ノ ハ エノキ も エノキ と 同様 に 人 為 的 な 影響 の 大 き 
な 地域 を 好む 結果 が 得 ら れ , 福田 ・ 中 内 (2010) (CLAM 
告 と 一 致し た . 一 方 で , エゾ エノキ は 気候 が 重要 な 要因 と 
し て 選択 され (#1), 生息 適地 も エノキ と は 異な る 分 布 を 
示し た . た だ し , AUC の 値 も 低かっ た こと か ら . モデ ル の 
信頼 性 は あま り 高 く な いと 判断 され る . 地方 版 レ ッ ド デー 
タブ ッ ク に よる と エゾ エノキ は 大 阪 府 で は 絶滅 (大 阪 府 . 
2000), 北海 道 と 鹿児島 県 で は 人 準 絶滅 危 恨 種 に 指定 され て 
BY (北海 道 環境 生活 部 , 2001 : 鹿児島 県 環境 生活 部 環境 
保護 課 、 2003)、 本 研究 に お いて も 全体 的 に 分 布 確率 が 低く 
推定 され て いる こと か ら (M4), 今回 の デー タ あ る い は 手 
法 で は 精度 の 高い モデ ル の 構築 が 難し い の か も し れ な い . 


本 研究 に よっ て 日 本 に お ける 外来 アカ ボ シ ゴ マダ ラム の 分 布 
は 気候 に よる 制限 を 受け に くい こと が 示唆 され , 定着 に 
は 寄主 植物 の 分 布 が 重要 な 役割 を 果たし て いる と 考え られ 
る . 寄主 植物 の 解析 結果 か ら 判断 する と , 外来 アカ ボ シ ゴ 
マダ ラ の 潜在 的 な 生息 適地 は , 現在 の 定着 確認 地点 より は 
る か に 広い こと が 示唆 され る . すなわち , PRC PUR IC Eo 
て 分 布 が 拡大 する 余地 が 多く 残さ れ て いる と いう こと で あ 
る . 本 研究 の 解析 に 用 いた エノキ 類 の 分 布 情報 は 植生 調査 
が 基 に な っ て いる . その た め , 公園 や 街路 樹 , 庭木 と し て 


植わっ て いる 植 栽 由来 の 個体 の 分 布 は 把握 で き て いな い . 
さら に , オオ ムラ サキ の 保全 活動 の 一 環 と し て . エノキ や 
エゾ エノキ が 植 教 され て いる 事例 も ある (高橋 , 1996 : 小 
林 , 2003). その た め , 実際 の エノキ 属 植 物 の 分 布 は 本 研究 
に よる 予測 より 広い 可能 性 が ある . この よう な 植 栽 は , 緑 
化 や 保全 の 観点 か ら 大 き な 意 義 が ある 一 方 で , 外来 生物 で 
ある アカ ボ シ ゴ マダ ラ の 分 布 拡大 を 意図 せ ず 促進 し て し ま 
4 可能 性 も ある . 


外来 アカ ボ シ ゴ マダ ラ の 在 来生 態 素 へ の 影響 と し て は , ほ 
ぼ 同 じ ニ ッ チ に 生息 する オオ ムラ サキ や ゴマ ダラ チョ ウ と 
の 競合 (松井 .2010 : 矢 後 .2012) や 交雑 ・ 繁 殖 阻害 (石垣 . 
2009 : KA, 2011 : 高 柔 .2012) が 懸念 され て お り , 本 研究 
の 結果 は 日 本 の 広い 範囲 が この よう な 危険 性 を 有 し て いる 
こと を 示唆 し て いる . さら に , 南西 諸島 に は クワ ノ ハ エ ノ 
キ や エノキ の 生息 適地 が 多く 存在 し て いた ( 図 4). 奄美 諸 
島 に は 固有 亜 種 で ある 在 来 ア カボ シゴ マダ ラ が 生息 し て お 
V, も し 外来 アカ ボ シ ゴ マダ ラ が 侵入 し た 場合 , 交雑 に よっ 
て 遺伝 的 汚染 が 発生 し て し まう 危険 性 が ある . 放 蝶 に よる 
新た な 分 布 拡大 を 阻止 すみ こと は も ちろ ん の こと , すでに 
定着 し て いる 外来 アカ ボ シ ゴ マダ ラ 集 団 や 寄主 植物 を どの 
よう に 管理 し て いく の か 十分 に 検討 し て いく 必要 が ある だ 
AY. 同時 に , 本 研究 に よる 外来 アカ ボ シ ゴ マダ ラ の 潜在 
的 生息 適地 の 広 さ を 考え る と , AM (2010) や 宇 式 (2010, 
2011, 2012) で 行わ れ て いる よう な 在 来生 態 系 へ の 影響 の 
調査 研究 を さら に 進め る 必要 が ある . 


放 蝶 ある い は 非 意図 的 な 導入 に よる 新規 の 定着 を 避け る 
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表 1. 


エノキ 属 3 種 に お ける ロジ ステ ィ ッ ク 回 帰 分 析 の 結果 . 
ば れ た 最適 モデ ル に 含ま れる 変数 の み 回 帰 係数 を 記し た . 


標準 化 回 帰 係数 標準 誤差 


モデ ル 選 択 に よっ て 選 





説明 変数 





エノキ エゾ エノキ クリ ワリ ノ ハ エノキ 
(定数 項 ) -2.863 + 0.044 -5.160 エ 0.144 -10.368 エ 0.897 
温 量化 指数 —1.238 + 0.359 
最低 気温 0.196 + 0.036 1.270 + 0.430 2.021 + 0.246 
夏 李 降水 量 ~0.070 + 0.048 
冬 李 降水 量 -0.101 + 0.044 0.268 + 0.083 
地形 的 湿潤 指数 -0.238 エ 0.161 -0.598 + 0.408 
丘陵 地 -0.139 +0.036 
山地 0.428 + 0.250 
景観 不 均一 性 0.148 + 0.075 0.696 + 0.388 
農地 率 0.527 + 0.059 0.907 + 0.381 
都市 率 0.428 + 0.041 1.136 + 0.254 


エノキ 





ai 生息 適合 度 
MRE 0.2262 - 0.4503 
デー Pe 0.1136 - 0.2262 
0.057 - 0.1136 
0 - 0.057 


生息 適地 









クリ ノ ハ エノキ 


生息 適合 度 
WE 0.2014 - 0.8582 
EEN 0.0472 - 0.2014 
0.0111 - 0.0472 
0- 0.0111 


生息 適地 








生息 適合 度 
ME 0.0314 - 0.0676 


Fea 0.0146 - 0.0314 
0.0068 - 0.0146 
0 - 0.0068 


生 








図 4. 最適 モデ ル に 基づく エノキ 属 3 種 の 生息 適合 度 . 凡例 は 分 布 確率 を 直し て お り , 


周 値 以上 の 適合 度 を 有する 地域 は . COR 








の 生息 に 適し て いる と 判断 され る . 
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よう に 促す た め に も ゃ , 侵入 前 ある い は 侵入 初期 の 段階 で 
分 布 解析 を 実施 し , 外来 召 の 管理 に 応用 する こと が 重要 
で あろ う . 本 研究 の 結果 を 見 る と , 外来 アカ ボ シ ゴ マダ 
ラ の 現状 も 事前 に 十分 に 予期 で きた と 考え られ る . AF 
完 の よう な アプ ロー チ は , 原産 地 の 標 本 情報 が あれ ば 容 
易 に 実行 する こと が で きる . 条件 に も よる が , Maxent の 
場合 は 数 十 地 点 の 標本 情報 で も 比較 的 頑健 な 結果 が 得 ら 
れる こと が 示唆 され て いる (Hernandez et al., 2006 ; Wisz 
et al., 2008). また , 行政 な ど が 組織 的 に 広域 で 収集 し て 
いる 生物 の 分 布 情報 を 使用 する こと が で きれ ば , FER 
物 の 視点 か ら も 解析 を する こと が で きる . た だ し , 現状 
で は 標本 情報 や 分 布 情報 の 公開 や 共有 は 発展 途上 に ある . 
曲 本 情報 な どの 分 類 学 研究 の 成果 は 外来 種 問題 に 有用 な 
情報 で ある こと は 疑い な く , GBIF (Global Biodiversity 
Information Facility : 地球 規模 生物 多様 性 情報 機構 ) な ど 
へ の 登録 も る 進め られ て いる が , 十分 と は 言い 難い . また , 
公園 や 道路 脇 な ど に 植 裁 され た 植物 の 情報 は , それ ら を 
寄主 と する 昆虫 の 外来 種 問題 を 考え る 際 に 重要 で ある が , 
情報 が 公開 され て いな いこ と も 多い . その た め , これ ら 
の 情報 が シー ムレ ス に 利用 で きる 形 で 共有 され る 枠組 み 
が 作ら れる こと が 望ま し い . これ まで に 数 多く 蓄積 され 
て きた 貴重 な 標本 情報 や 分 布 情報 が 多く の 人 々 に 共有 さ 
れ , 外来 生物 対策 の よう な 応用 的 な 局面 に 広く 活用 され 
る こと が 期待 され る . 


H 辞 

環境 省 生 物 多様 性 セン ター に は , 第 6・7 回 自然 環境 保全 
基礎 調査 に お ける 植生 調査 デー タ の 使用 許可 を いた だ い 
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Summary 


Despite the existence of native populations of Hestina assimilis 
(Linnaeus) and its allied species in Japan, a continental population 
of H. assimilis has recently been introduced into Kanagawa Pref. 
and is spreading its distribution range, mainly in the Kanto region. 
This introduced population seems to have a negative effect on the 
native ecosystem. To consider the appropriate management of 
range-expanding introduced butterflies, it is important to project the 
potential distribution by a species distribution model. In this study, 
we estimated the potential habitat of the introduced population of H. 
assimilis in Japan, based on climate and host-plant information. 


From specimen records in China and South Korea, we obtained 49 
coordinates of H. a. assimilis in its native range. We constructed 
a Maxent model which explains the presence of H. a. assimilis by 
climate factors (i.e. minimum temperature of the coldest month 
and minimum monthly mean of relative humidity in winter season) 
in China and South Korea. The climatic potential habitat of H. a. 
assimilis in Japan was projected by extrapolating this model to 
Japan. We estimated the potential distribution of three Celtis plants, 
which are host plants of H. a. assimilis, using the result of vegetation 
surveys over a wide area of Japan. To obtain the potential habitat of 
Celtis, logistic regression analysis was applied to model the presence 
or absence of Celtis based on climate (i.e. warmth index, monthly 
mean of the daily minimum temperature of the coldest month, and 
total summer and winter precipitation), topographic factors (i.e. 
topographic wetness index, the presence or absence of hills, and the 


presence or absence of mountains) and land use (i.e. the proportion 
of agricultural area and urban area and landscape heterogeneity). 


On the basis of the Maxent model, many areas in Japan had a 
suitable climate for H. assimilis. Thus, in colonization of H. a. 
assimilis in Japan, the distribution of host plants (i.e. Celtis) has 
played the critical role. Celtis sinensis Pers. and Celtis boninensis 
Koidz. preferred artificial environments, while C. sinensis had a 
wide range of suitable habitat. No reliable model was obtained for 
Celtis jessoensis Koidz. Our results indicate that the potential habitat 
of H. assimilis in Japan is much wider than the present distribution 
area. This indicates that the potential impact on the native 
ecosystem is also much broader. Because there may be competition, 
hybridization, and interspecific promiscuity with its allied species, 
Hestina persimilis (Westwood) and Sasakia charonda (Hewitson), 
which inhabit almost the same niche, as a result of range expansion 
of H. a. assimilis, further research is required on the influence on 
the native ecosystem. Our approach is very useful and effective 
where specimen records in the native range are available. Moreover, 
potential habitat analysis in terms of host plant could also be 
conducted if sufficient distribution information on plants collected 
in the wide area is available. It is expected that the publication 
and sharing of distribution information on organisms will promote 
appropriate future wildlife management, such as the management of 
introduced butterflies. 
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